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Abstract 

The experiments of Farmer, Hanson, and Wynne^ 
showed that use of a supercritical airfoil can 
adversely affect wing flutter speeds in the tran- 
sonic range. Inasmuch as adequate theories for 
three-dimensional unsteady transonic flow are not 
yet available, the modified strip analysis published 
by Yates in 19583 has been used to predict the 
transonic flutter boundary for the supercritical 
wing tested by Fanner, Hanson, and Wynne. The 
steady-state spanwise distributions of section 
lift-curve slope and aerodynamic center, required 
as input for the flutter calculations, were obtained 
from pressure distributions measured by Harris. 2 
The calculated flutter boundary is in excellent 
agreement with experiment in the subsonic range. 

In the transonic range, a "transonic bucket" is 
calculated which closely resembles the experimental 
one with regard to both shape and depth, but it 
occurs at about 0.04 Mach number lower than the 
experimental one. 


Nomenclature 

a nondimensional distance from midchord to 

’ section aerodynamic center measured 

perpendicular to elastic axis, positive 
rearward, fraction of semi chord b 
b semichord of wing measured perpendicular 

to elastic axis at spanwise reference 
station n = 0.75 (bp = 0.14948 m) 

C. section lift-curve slope for a section 

a,n perpendicular to elastic axis 

Cp pressure coefficient 

H translational displacement of wing at 

elastic axis, positive downward, 
fraction of reference semi chord bp 
k reduced frequency based on spanwise 

reference station (n = 0.75) and on 
velocity component normal to elastic 
axis, bpO)/V cos Aea 

M freestream Mach number 

mp mass of wing per unit span at spanwise 

reference station (n = 0.75) 

(mp = 2.57 kg/m) 

q freestream dynamic pressure 

V freestream speed 

X streamwise distance from local leading 

edge, fraction of local chord 

Y spanwise station, fraction of semi span 
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a angle of attack 

n nondimensional coordinate measured from 

wing root along elastic axis, fraction 
of elastic axis length 

0 torsional displacement of wing about 

elastic axis, positive leading edge up 
Aea sweep angle of elastic axis (Aea = 40.03®) 

Up mass ratio based on spanwise reference 

station (n = 0.75), mp/7Tpbp2 

p freestream density 

oj circular frequency of vibration 

u)p reference frequency, frequency of first 

uncoupled torsional mode of wing 
(o)p = 227.23 rad/sec) 


Introduction 

The experiments of Farmer, Hanson, and Wynne, 1 
which determined the flutter characteristics of a 
supercritical wing and a dynamically similar con- 
ventional wing having the same planform, showed 
that use of a supercritical airfoil can adversely 
affect wing flutter speeds in the transonic range. 
Although it was evident from aerodynamic data 
available more than ten years ago2 that this adverse 
effect should occur, the results of Ref. 1 aroused 
considerable concern, which in turn led to renewed 
interest in predicting the flutter characteristics 
of supercritical wings. Inasmuch as conventional 
lifting-surface theories do not lead to accurately 
predicted transonic flutter characteristics, and 
because adequate theories for three-dimensional 
unste^idy transonic flow have not yet been developed, 
the modified strip analysis, first published by 
Yate$3 in 1958, has been used here to predict the 
transonic flutter boundary for the supercritical 
wing (Fig. 1] which was investigated in Ref. 1. 
McGrew et an have also calculated flutter charac- 
teristics for this wing but by use of subsonic 
doublet-lattice aerodynamic in combination with an 
extensive set of "weighting factors." 

The modified strip analysis used here has 
given good flutter results for a broad range of 
swept and unswept wings at speeds up to hypersonic, 5 
including effects of wing thicknessS»7 and angle of 
attack. 8 in particular, it was used successfully 
to calculate transonic flutter characteristics for 
some swept wings with conventional airfoils.^ 

In addition to results calculated for com- 
parison with the experimental data of Ref. 1, 
results are also shown herein which more broadly 
define the flutter characteristics of this wing in 
terms of the variations of flutter parameters with 
mass ratio at constant Mach number and with Mach 
number at constant mass ratio. 
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Flutter Model Parameters 

The geometric, elastic, and inertial properties 
used in the present calculations are measured 
values for the supercritical -wing flutter model of 
Ref. 1. The model geometry is shown in Fig. la. 

The frequencies and mode lines of the first six 
measured natural vibration modes are presented in 
Fig. lb, and the modal deflections are given in 
Fig. 2. The corresponding generalized masses were 
determined by the method of displaced frequencies.lO 

These six measured modes were used in all of 
the flutter calculations shown herein. Some col- 
lateral flutter calculations, made with up to ^ 
twelve vibration modes calculated with the NASTRAN* 
finite-element structural analysis, indicated that 
six modes were sufficient to converge flutter 
speeds and frequencies within 2 percent. Flutter 
calculations with the measured modes, however, 
appeared to converge more rapidly than those with 
NASTRAN modes. Consequently, the six-measured -mode 
results presented here are considered to be con- 
verged within about 1 percent. 


Flutter Analysis 

Method 

The modified strip analysis^ used here is 
formulated for wing strips oriented normal to the 
elastic axis and is based on stripwise application 
of Theodorsen-type aerodynamics^^ in which the 
lift-curve slope of 2 t 7 and aerodynamic center at 
quarter chord are replaced, respectively, by the 
lift-curve slope and aerodynamic center for the 
same strip of the three-dimensional wing at the 
appropriate Mach number. The downwash collation 
point, where the downwash induced by the aerodynamic 
load is set equal to the kinetmatic downwash, is 
modified accordingly. The circulation function is 
modified for compressibility by use of two-dimen- 
sional unsteady compressible-flow theory. 12 
Further description and discussion of this method 
are contained in Refs. 3, 5, 9, and 13. 

Aerodynamic Parameters 

The required spanwise distributions of section 
lift-curve slope and aerodynamic center were 
obtained from steady-state surface pressure measure- 
ments made by Harris^ in the Langley 8-Foot Tran- 
sonic Pressure Tunnel and in the Langley 16-Foot 
Transonic Tunnel. Upper-and lower-surface pressures 
were measured along streamwise wing sections at 
six stations along the semispan (Figs. 3 and 4). 
Since the flutter tests were conducted at angles 
of attack near zero, the pressure data for the two 
angles of attack nearest zero were selected for 
use in the flutter calculations. The boundary- 
layer transition strips on the pressure model were 
represented also on the flutter model. 

Values of section lift-curve slope and section 
aerodynamic center for 40 wing sections equispaced 
along the elastic axis and oriented perpendicular 
to the elastic axis were required as input for the 
flutter calculations. These values were obtained 
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by differencing the corresponding values of section 
lift coefficient and pitching-moment coefficient 
for the two angles of attack shown on each part of 
Figs. 3 and 4, The section lift and moment coef- 
ficients were calculated by 8-point Gaussian 
integration of the lifting pressure along each of 
the 40 wing sections 

The procedure for obtaining values of lifting 
pressure at the eight quadrature abscissas on each 
of the 40 wing sections from the measured pressures 
is described in the following discussion and sketch. 



® Two of six streamwise lines along which 
pressures were measured 

(D One of eight points on each streamwise line at 
which pressures were interpolated 

(3) One of 40 wing sections on which lift and 
moment coefficients were required 

(4) One of eight lines of constant quadrature 
abscissa 

(5) One of eight quadrature points on each of 
40 wing sections 


Since measured pressures were available only along 
six streamwise lines (see sketch), it was necessary 
to interpolate the upper- and lower-surface pres- 
sures streamwise (Figs. 3 and 4) for the eight 
Gaussian-quadrature abscissas using the spline 
curves shown in these figures. For each of the 
quadrature abscissas, the resulting six values of 
lifting pressure (lower-surface pressure minus 
upper-surface pressure) were spline-interpolated 
from root to tip, and values were thus obtained 
at the required 40 wing sections. 

The distributions of section lift-curve slope 
and section aerodynamic center thus obtained are 
shown in Figs. 5 and 6. The change in the slopes 
of the curves inboard is caused by use of the 
dashed planform in Fig. 1 as reference while pre- 
serving the section lift and moment values for the 
wing with glove. The spline-curve extrapolations 
beyond the aft-most experimental points in Figs. 3 
and 4 generally indicate nonvanishing lifting 
pressures at the wing trailing edge. To examine 
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the extent to which the calculated flutter charac- 
teristics were affected by the pressure-data fair- 
ing near the trailing edge, some collateral cal- 
culations were made with a Kutta condition imposed 
to force the lifting pressure to zero at the 
trailing edge. The resulting changes in calculated 
flutter characteristics were negligible. 

Pressure data^ from the tests in the 8-Foot 
Transonic Pressure Tunnel for Mach numbers up to 
1.0 are available for two levels of dynamic pres- 
sure at each Mach number. Consequently, it is 
possible to assess the effect on aerodynamic param- 
eters (and hence on flutter) of static aeroelastic 
deformation of the pressure model, together with an 
accompanying change in Reynolds number.* The 
comparison in Fig. 5f is typical. Increasing the 
dynamic pressure reduces the lift-curve slope and 
moves the aerodynamic center forward, especially 
over the outboard portions of the wing. These 
changes have opposing effects on flutter speed, 
and the resulting changes in calculated flutter 
characteristics are quite minor. Nevertheless, 
the lower value of dynamic pressure was selected 
for all of the flutter calculations presented 
herein because those pressure data should be less 
affected by deformation than those obtained at the 
higher dynamic pressure. The effect of the asso- 
ciated change in Reynolds number is considered to 
be quite small since the previously indicated 
changes in aerodynamic parameters caused by chang- 
ing dynamic pressure are in a direction opposite 
to that anticipated from Reynolds-number change 
alone. 

Mass Ratio 

The mass ratios used in flutter calculations 
made for direct comparison with the experimental 
flutter data of Ref. 1 were taken from a curve 
faired through the experimental mass ratios (Fig. 7). 
Inasmuch as no flutter data were obtained at Mach 
numbers higher than 1.055, the mass ratio for that 
Mach number was used for all higher Mach numbers. 

In addition, mass ratio was varied parametrically 
up to a value of 98 for each Mach number in order 
to trace out "slices" of a flutter surface (flutter- 
speed index as a function of Mach number and mass 
ratio) such as that discussed in Appendix C of 
Ref. 6. 


Results and Discussion 

Values of flutter-speed index arid flutter- 
frequency ratio calculated with mass ratios taken 
from the curve in Fig. 7 are compared with the 
experimental data^ in Figs. 8 and 9, respectively. 
Agreement between calculated and measured flutter 
boundaries (Fig. 8) is excellent in the subsonic 
range. In the transonic range, a "transonic bucket" 
is calculated which closely resembles the experi- 
mental one with regard to both shape and depth. 
However, the calculated bucket occurs at about 


*Reynolds numl)ers of 2.2 x iq 6 and 3.3 x iq 6 at 
M = 0.90 (based on mean geometric chord) are 
representative for the 8-foot tunnel tests. The 
corresponding Reynolds number for the flutter 
test was near the higher value. 


0.04 Mach number lower than the experimental one. 

The reason for this difference is not known with 
certainty. There is some indication, however, 
that the difference may be associated with flutter- 
model scale effect in the wind tunnel. A limited 
number of flutter points were measured with a 
second model of the wing with conventional airfoil* 
but 0.4 as large as that reported in Ref. 1 (and 
0.4 as large as the supercritical-wing model). The 
data for the smaller model showed the bottom of 
the transonic bucket to be shifted down about 0.04 
in Mach number relative to that for the larger 
model. It is not likely that this difference was 
caused by the difference in blockage because 
blockage would be expected to cause the larger 
model to encounter transonic effects at lower 
Mach numbers than the smaller model. It is noted 
that, relative to tunnel dimensions, the pressure 
model is intermediate in size between the two con- 
ventional-airfoil flutter models and smaller than 
the supercritical wing model. 

The calculated flutter frequencies (Fig. 9) 
are slightly higher than experimental values up to 
Mach number 0.95, as was the case in Ref. 4. The 
reason for the somewhat erratic variation between 
Mach numbers 0.95 and 1.00 (square symbols in 
Fig. 9) is not clear. In that Mach number range, 
however, aerodynamic data are available from both 
the 8-foot tunnel and the 16-foot tunnel. Although 
the calculated flutter speeds are in good agreement 
with each other at the two duplicated Mach numbers 
(M = 0.95 and M = 0.99) (Fig. 8), the flutter fre- 
quencies obtained from the 16-foot tunnel data 
(diamond symbols) are higher and farther from 
flutter experiment than those calculated with 8-foot 
tunnel data (square symbols). The cause of this 
difference is entirely aerodynamic, and this may 
indicate that there is also a model scale effect on 
flutter frequency. Relative to tunnel size, the 
pressure model in the 16-foot tunnel was smallest; 
the same pressure model in the 8-foot tunnel was 
intermediate; and the flutter model in the Tran- 
sonic Dynamics Tunnel was largest. Thus, if 
extrapolation of flutter frequency with respect to 
model size is valid, lower calculated flutter 
frequencies would be expected if the pressure 
model were the same size as the flutter model rela- 
tive to tunnel dimensions. This is, of course, 
speculative and empirical since no physical 
mechanism to produce this effect is postulated. 

The variations of flutter-speed index, flutter- 
frequency-ratio, and reduced frequency with mass 
ratio (Fig, 10) are quite conventional and confirm 
that the upturn in the flutter speed and flutter 
frequency (Figs. 8 and 9) as Mach number decreases 
to 0.25 is caused by the accompanying decrease in 
mass ratio. This effect is illustrated more 
clearly in Fig. 11 where the calculated flutter 
characteristics from Figs. 8 and 9 are compared 
with corresponding values calculated for a constant 
mass ratio of 27.41 which is the value for Mach 
number 0.80 from Fig. 7. 


Concluding Remarks 

Flutter calculations have been made by modified 
strip analysis for a supercritical wing model for 
which experimental flutter data and steady-state 
pressure distributions were previously available. 
Use of these pressure data to generate aerodynamic 
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input for the flutter calculations produced a 
flutter boundary that is in excellent agreement 
with experiment in the subsonic range. In the 
transonic range, a "transonic bucket" was calcu- 
lated which closely resembles the experimental one 
with regard to both shape and depth, but it occurs 
at about 0.04 Mach number lower than the experi- 
mental one. Some evidence indicates that this 
shift may be related to differences in model size 
relative to tunnel size for the pressure model 
and the flutter model, but this is not conclusive. 

Nevertheless, the good results herein for the 
supercritical wing and good results obtained pre- 
viously for swept wings with conventional airfoils 
at subsonic, transonic, and supersonic speeds 
indicate that the modified strip analysis in 
conjunction with aerodynamic parameters from 
steady-state experiments, is still useful for 
transonic speeds in the absence of validated non- 
linear, three-dimensional, unsteady aerodynamic 
analysis methods. 
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MODE 1, 33.24 rad/ sec MODE 2, 113.6 rad/ sec 





(b) Node lines of measured vibration modes. 
Fig. 1 Concluded. 
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M = .50, 


Fig. 3 


O UPPER SURFACE 
□ LOWER SURFACE 
O a = -.08 deg. 
Cf a = 2.02 deg. 



M = .80. q = 31.1 kPa 

Fig. 3 Continued. 
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Fig. 3 


O UPPER SURFRCE 
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M = .95, 


Fig. 3 


o UPPER SURFACE 
□ LOWER SURFACE 
O « — --07 deg. 
CT a = 2.32 deg. 








0 


.2 


.4 


.6 


.8 1.0 


X 


(h) M 

Fig. 


1 . 00 , 

Concli 


14 


O UPPER SURFACE 
□ LOWER SURFACE 
O « — .05 deg. 

O' a = 2.39 deg. 
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q = 28.7 kPa 
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(a) M = .95. q = 35.4 kPa 


Fig. 4 Pressure distributions measured in 16-Foot Transonic Tunnel 



M = .99. q = 36.7 kPa 

Fig. 4 Continued. 




Fig. 4 Continued. 
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Fig. 4 Concluded. 
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(e) M = .95, q = 28.7 kPa 
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(f) M = .98 

Fig. 5 Continued. 
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(a) M = .95, q = 35.4 kPa 
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(b) M = .99. q = 36.7 kPa 


Fig, 6 Aerodynamic parameters for flutter analysis from 16- foot tunnel data. 
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(c) M =1.02, q = 37.6 kPa 



(d) M =1.20, q = 41.6 kPa 


Fig. 6 Concluded. 
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Fig. 7 Mass ratios for flutter experiments. 



Fig. 8 Comparison of calculated and measured flutter speeds. 
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SOLID SYMBOLS INDICATE LOW DAMPING 
O — FLUHER EXPERIMENT (REF. 1) 

□—AERO. PARAMETERS FROM 8-FOOT TUNNEL 
O— AERO. PARAMETERS FROM 16-FOOT TUNNEL 
KERNEL FUNCTION (REF. 1) 


Fig. 9 Comparison of calculated and measured flutter frequencies. 
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(a) FLUTTER SPEED INDEX 

Fig. 10 Variation of calculated flutter characteristics with mass ratio. 
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Fig. 11 Comparison of calculated flutter characteristics for constant and varying mass ratio. 
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